arXivrl 506.03882V 1 [astro-ph.HE] 12Jun2015 


Mon. Not. R. Astron. Soc. 000,^]^ (2015) Printed 15 June 2015 


(MN style file v2.2) 


Revisiting the Dynamical Case for a Massive Black Hole in 
ICIO X-1 

Silas G. T. Laycock^*, Thomas J. Maccarone^ and Dimitris M. Christodoulou^ 

^Department of Physies, University of Massaehusetts Lowell, MA, 01854, USA 

‘^Department of Physies, Texas Teeh University, Lubboek TX, 79409, USA 

^Department of Mathematieal Seienees, University of Massaehusetts Lowell, MA, 01854, USA 


15 June 2015 


ABSTRACT 

The relative phasing of the X-ray eclipse ephemeris and optical radial velocity (RV) 
curve for the X-ray binary ICIO X-1 suggests the He[A4686] emission-line originates 
in a shadowed sector of the stellar wind that avoids ionization by X-rays from the 
compact object. The line attains maximum blueshift when the wind is directly toward 
us at mid X-ray eclipse, as is also seen in Cygnus X-3. If the RV curve is unrelated to 
stellar motion, evidence for a massive black hole evaporates because the mass function 
of the binary is unknown. The reported X-ray luminosity, spectrum, slow QPO, and 
broad eclipses caused by absorption/scattering in the WR wind are all consistent with 
either a low-stellar-mass BH or a NS. For a NS, the centre of mass lies inside the WR 
envelope whose motion is then far below the observed 370 km/s RV amplitude, while 
the velocity of the compact object is as high as 600 km/s. The resulting 0.4% doppler 
variation of X-ray spectral lines could be confirmed by missions in development. These 
arguments also apply to other putative BH binaries whose RV and eclipse curves are 
not yet phase-connected. Theories of BH formation and predicted rates of gravitational 
wave sources may need revision. 

Key words: accretion, accretion disks—circumstellar matter—stars: neutron—X- 
rays: binaries—black holes 


1 INTRODUCTION 


1.1 High-Mass X-Ray Binaries 


ICIO X-1 is a high-mass X-ray binary (HMXB) containing a 

l38__ 


wind of its Wolf-Rayet (WR) companion ( [Prestwich et al. 
[2007 ). The binary period is 35 hrs and deep eclipses are seen 
in the X-ray light curve. This Letter addresses the strong 
doubts cast by recent observations on the massive black hole 
(BH) interpretation for this system. A substantial downward 
revision of the mass of the compact object will impact re¬ 
search on the evolution of massive stars and BHs of all sizes 


(Barnard et al.||2008 

Belczynski et al.||2010 Carpano et al. 

2008 Park V Riccoti 

2011), long gamma-ray bursts and BH- 


leTaLllMTI l, as well as the so-called Ultra-Luminous X-ray 


(ULX) sources ( Soria||2QQ7 Gladstone et al. 2009 Feng V 


Soria|[20TT Medvedev V Putanen||20T3 ). In particular, the 


noted preference for massive BH candidates to occur in low 
metallicity environments would lose its central exemplar. 


E-mail: silas_laycock@uml.edu 


Binary systems containing X-ray emitting compact objects 
accreting matter from their companion stars are thought to 
host both neutron stars (NSs) and black holes. The compact 
objects themselves are invisible due to their very small sizes 
and correspondingly low intrinsic luminosities. Their exis¬ 
tence is inferred from the properties of the X-ray emission 
and from the doppler-derived motions of the companions as 
the stars orbit around their centers of mass. Spectral lines 
originating at the surface of the primary are observed re¬ 
peatedly over the course of one or more orbital cycles to 
establish the amplitude of the radial velocity (RV) curve. 
Tidal distortion of the primary (in close systems) also leads 
to a periodic optical brightness modulation, known as ellip¬ 
soidal variation, which can be seen in the photometric light 
curve. Modeling of these variations can yield an independent 
measurement of the orbital period and inclination of the sys¬ 
tem. Since it is only the primary that is directly observed, 
the mass function is combined with independent estimates 
of its mass to constrain the mass of the unseen companion. 
These techniques have yielded a growing list (e.g., Casares 
V Jonker|2014| Remillard V McClintock 2006) of binary sys¬ 


tems where the compact object mass is significantly above 
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the theoretical limit of 2.9 M© for neutron stars. Such sys¬ 
tems are generally (although not universally) accepted to 
contain BHs. 

Another class of objects thought to contain stellar-mass 
BHs are the ULX sources (Liu et al. 2013 Motch et al.|2014 ). 
Their X-ray luminosities (10^^ — 10"^^ erg s“^) exceed the Ed- 
dington limit for a NS by 1 — 3 orders of magnitude. Most 
ULXs are too distant for RV or photometric studies by ex¬ 
isting instruments to reveal their dynamical parameters, so 
the argument for stellar-mass BHs in these systems derives 
mainly from their X-ray properties. 

The X-ray binary ICIO X-1 occupies an important niche 
as a link between BH-HMXBs and ULXs due to its inter¬ 
mediate luminosity (~10^®erg s“^) and the fact that, at 
a distance of 660 kpc, it is one of the few extragalactic 
systems close enough to be readily observable by current 
generation telescopes. The landmark work of [Silverman 


Filippenko (2008) revealed that ICIO X-1 appears to have 


an RV-derived mass function that requires a BH of 23-32 
M©. SF noted that their RV curve had not been phase- 
aligned with the X-ray eclipses since the X-ray period was 


not then known to sufficient precision. Others (e.g. Tutukov 
et al.|20l'3 ) have pointed out that this caveat also applies to 


NGC300 X-1 which is in many regards a twin to ICIO X-1 
in terms of its WR companion, similar orbital period, and 
X-ray luminosity. 

1.2 Wind Ionization and Shadowing 

In the most massive and hottest stars, the stellar photo¬ 
spheres are highly ionized so that no suitable absorption 
lines can be found and emission lines are used instead. 
The most suitable ion species is He H because it can ex¬ 
ist at the temperatures found in the photospheres of O stars 
(> 3 X lO^K), where it gives rise to the Pickering series of 
spectral lines. Emission lines arise in low-density gas, so the 
region probed by the He H lines can in principle extend to 
many stellar radii above the surface. The density of stellar 
winds falls with distance r approximately as 1/r^, so most 
of the equivalent width of the line should originate at lower 
altitudes. Provided that the wind is isotropic (or at least 
azimuthally symmetric), the He H line should still trace the 
orbital motion of the star. 

At the very high X-ray fluxes seen in high luminosity 
XRBs, the wind will be completely ionized out to a distance 
far exceeding the orbital separation of the binary. This fact 
has not been generally recognized as affecting the mass de¬ 
terminations obtained via the RV-curve method. IMaccaronel 


et al. (2014) have recently pointed out that a model intro¬ 


duced by van Kerkwijk et al. (1996) for Cyg X-3 might ap¬ 


ply more widely. In this model. X-ray photoionization of 
the wind suppresses line emission throughout the binary or¬ 
bit. Only in a small sector of the wind shielded by the pri¬ 
mary will the gas be able to maintain a “normal” ionization 
state. In ICIO X-1, the ionization parameter has a large value 


{Lxlndr ~ 10^ ) at the location of the WR star (Barnard et 


[ar]|2014| . The wavelength and shape of the observed emis¬ 
sion lines will be a convolution of the wind velocity Vwind, 
the orbital motion vorbits and the stellar rotation Vrot- 

In the limit that Vwind » Vorbit,Vrot, the observed RV 
curve will have an amplitude of Vwind sinz, where i is the in¬ 
clination of the orbital plane to the line of sight. The phase 


of the resulting RV curve will be as follows: at mid X-ray 
eclipse, when the WR star lies between us and the X-ray 
source, the shielded wind sector is pointing toward the ob¬ 
server, as illustrated in Figure The wind velocity vector 
then lies along our line of sight, producing the maximum 
blueshift of the emission line. At quadrature, when the WR 
star’s projected radial velocity is greatest, the wind in the 
shielded sector is orthogonal to the line of sight, producing 
zero doppler shift. Finally at superior conjunction, the emit¬ 
ting portion of the wind is expanding away from us, leading 


to the maximum redshift. Hanson et al. (2000) were able to 


map this cool structure in Cyg X-3 from IR spectra using 
doppler tomography. 

The actual phase relationship between X-ray eclipses in 
ICIO X-I and the RV curve of the optical counterpart has 
been recently explored by Laycock et al. (2015). After con¬ 


structing a new ephemeris for the X-ray eclipses from a 10 
year series of Chandra and XMM-Newton observations, it 
was found that the published radial velocity data of [Silver^ 
man & Filippenko (2008) are offset by 1/4 of an orbit from 


the geometric alignment implied by the eclipses. The phase 
shift is such that the periodic doppler shift of the He [A4686] 
line does not pass though zero at mid-eclipse. Instead, it is 
the maximum blueshift that occurs at mid-eclipse which is 
expected to coincide with inferior conjunction of the WR 
primary. 

Variations in He H line-width are evident from Ta¬ 


ble 2 of Silverman V Filippenko (2008), which suggests the 


FWHM is largest when the RV is closest to zero, and small¬ 


est at maximum redshift. In van Kerkwijk et al. (1996) the 
line profile model consists of a weak broad contribution from 
the hot mostly ionized stellar wind, plus a narrow velocity- 
modulated line from the cool gas in the shadow which alter¬ 
nately passes in front of and behind the star, attaining its 
minimum width (and greatest strength) at maximum red¬ 
shift. Vwind in the line-forming region is substantially lower 
than the terminal wind velocity. jVilhu et al. (2009) used 
Chandra HETG spectra to map the ionization and velocity 
fields of the wind in Cyg X-3, largely confirming the pic¬ 


ture developed by van Kerkwijk et al. (1996) and Hanson et 


al. (2000) and showing that the wind is highly ionized out 


to many times the binary separation. The X-ray EeXXVI 
line RV curve appears to trace the compact object’s motion, 
with the required phase relationship to the IR emission and 
absorption features, resulting in a mass function consistent 


with that obtained by Hanson et al. (2000) 


The appeal of the van Kerkwijk et al. (1996) model is 


that it isolates a single phenomenon (photoionization) that 
seems likely to dominate interaction in IC 10 X-I; it is not 
unique in producing phase-shifts and newer models for Cyg 
X-3 successfully reproduce many elements of its more com¬ 
plex dynamics. The phenomena encompassed by those mod¬ 
els probably apply to a lesser extent in ICIO X-I since the 
orbital separation in Cyg X-3 is smaller, resulting in the 
compact object being engulfed by the stellar envelope. The 
flow of material in that system is influenced by both ion¬ 
ization and tidal interactions, leading to density structures 
such as plumes, wake-trails, and shocks; while ICIO X-I has 
so far been observed in a state analogous to the quiescent 
state of Cyg X-3. 
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1.3 Spectral Analysis of ICIO X-1 

ICIO is distant and so grating spectra have not yet been 
obtained. [Wang et al. (2005) fit the XMM-Newton EPIC 
(pulse height) spectrum with a multicolor disk blackbody 
model appropriate to a BH in the low hard state. Finding 
an inner accretion disk temperature of Tjn=l.lllo os 
at a radius of Rin—‘^^-Q km, disk inclination angle > 57°, 
subsolar metal abundance and Nh ~10^^cm“^. The point 
source is blended with an extended diffuse structure of effec¬ 
tive radius of 501^^ km, which appears to be an ionized X- 
ray nebula surrounding the compact object. Unfortunately, 
these early observations did not accumulate enough counts 
during eclipse to fully isolate the spectral components. By 
looking at the variation in hardness ratio (HR) during eclipse 
egress in the Chandra data, Laycock et al. (2015) found evi¬ 


dence of energy dependence and noted that the hardness ra¬ 
tio increases before the flux begins its climb at egress. Energy 
dependence is not expected in the eclipse profile of an X-ray 
point source passing behind a fully ionized absorber. The 
modulation would be dominated by free electron scattering 
which is energy independent over this wavelength range. At 
the more moderate ionization fraction expected deep in the 
core of the wind (and especially in the shadow of the WR 
star), single-electron ions of e.g. O, C, Si, and Fe will be 
present (as seen by |Vilhu et al. (2009) in Cyg X-3) and con¬ 
tribute significant soft X-ray absorption. 

A much longer XMM-Newton EPIC observation an¬ 


alyzed by Barnard et al. (2014) reveals evidence for an 


accretion-disk corona. Their Fig. 2 shows a clear HR mod¬ 
ulation during eclipse. The interpretation is that different 
parts of the emission region (disk and corona) become cov¬ 
ered by the absorber as the eclipse progresses. They avoid 
the word eclipse in favour of “dip” because the passage of the 
compact object through (and behind) the wind gives rise to 
varying amounts of scattering and absorption at all phases. 
Thus, what we call the eclipse is referred to by them as the 
“main dip.” The inner-edge temperature for their preferred 
accretion-disk spectral fit (Model V) is Tin = l.6S keV, and 
this component disappears during eclipse, while about 10% 
of the spatially-extended, comptonized corona component 
remains visible (“partial covering”). We note that Model V 


is not the best-fit model of Barnard et al. (2014). Model HI 


with Tin =0.2 keV and a compact corona was a better fit, 
but was rejected (in part) because it does not conform to 
the assumption of a (massive) stellar mass BH. The very 
large number of parameters required in these models results 
from the lack of information about the spectrum above 10 
keV, and the inability of current datasets to resolve the soft 
X-ray lines expected from hydrogen-like heavy ions in the 
absorbing core of the wind. 

In what follows, we adopt the model shown in Figure 
for optical line emission during X-ray eclipses and in § 2 we 
make the case for a low-mass compact object in ICIO X-1, 
possibly a NS although a low-mass BH cannot be ruled out. 
In § 3, we discuss the implications for other HMXBs and 
ULXs, and in § 4 we summarize our conclusions. 


2 A LOW-MASS COMPACT OBJECT 

If we in fact do not know the mass function of ICIO X- 
1, then the only hard evidence on which to base a model 


Table 1. Orbital Parameters for a NS or a BH in ICIO X-1 


Mwr 

Mx 

a 

Rwr 

vx 

R(s) 

Rwr 

Twr 

Mq 

Mq 

Rq 

km/s 

km/s 

Rq 

Rq 

hr 

35 

1.4 

18.99 

23.96 

599.10 

8.06 

2 

1.24 

35 

2.5 

19.18 

41.95 

587.32 

8.14 

2 

1.23 

17 

1.4 

15.13 

37.76 

458.57 

6.42 

1.5 

1.17 

17 

2.5 

15.42 

64.88 

441.15 

6.55 

1.5 

1.15 

35 

32 

23.27 

364.7 

398.9 

9.88 

2 

1.01 

17 

23 

19.59 

369.7 

273.3 

8.32 

1.5 

0.90 


Radius of shadowing body R(s) determined from observed eclipse 
duration Tmin=^ hr. Inferred duration for eclipse by the WR star 
is Twr- Table assumes inclination i = 90° 


for the system comes from the orbital period and the stel¬ 
lar parameters of the WR primary. The X-ray luminosity 
and eclipse duration are informative but they are subject 
to many more parameters (wind density, ionization, accre¬ 
tion efficiency, magnetic field), and consequently they can 
be used to support competing non-unique solutions. Here 
we analyze the plausibility of a much lower mass BH or a 
NS as the compact object. 

The mass of the WR star is reasonably well constrained 
by the optical classification of Clark X Crowther (2004) with 
a most probable value of Mwr^^^Mq and a range of 17-35 
Mq. For a neutron star companion, a commonly accepted 
range for the mass is Mns = 1.26 — 2.5 M©, where the lower 
bound is the Chandrasekhar limit for a white dwarf minus 
the binding energy lost in collapse to a NS, and the upper 
bound is the Tolman-Oppenheimer-Volkoff (TOV) limit for 
a NS. The highest reliable measurement is 2.0 Mq and cur¬ 
rent EOS allow 0.1-2.9 Mq. For BHs there are no physical 
bounds to the mass and any dynamically confirmed object 
significantly in excess of 2.5 Mq is regarded as a probable 
BH. 

Applying Kepler’s law with NS masses yields an or¬ 
bital separation of (ai + a 2 ) = IARq for the lowest mass 
combination (1.4 + I 7 M 0 ), and I 9 R 0 for the highest mass 
combination (2.5 + 35Mq). Naturally, the NS orbit shrinks 
now to 14-19^0 relative to the 19-23^© range used in prior 
investigations that assumed a 23-32 M© BH (see Table 0. 
At these closer separations, the stellar wind density encoun¬ 
tered by the compact object increases and the primary ap¬ 
proaches filling its Roche lobe. The size of the Roche lobe 
of the primary for all of these combinations of WR star and 
NS masses is ri 0.6(ai + a 2 ) = S-11Rq. 

Analyses of the Chandra and XMM-Newton eclipse pro¬ 
files by Laycock et al. (2015) and [Barnard et "ar (2014), 
respectively, measured slightly different durations for mini¬ 
mum flux (5 hr vs. 5.2 hr), and ingress/egress phases (~1 
hr vs. ~2 hr). This appears to be due to the broader and 
more gradual eclipse seen by XMM-Newton. Differences in 
the method used to estimate the diameter of the WR star 
also led to small variations in the estimates of the size of the 
eclipsing body (4 vs. 5 for Rwr) . [Barnard et al. (2014) esti¬ 
mated the effective stellar radius Rwr seen by the X-rays to 
be the blackbody radius implied by the color temperature 
reported by Clark X Crowther (2004); while we used the 
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tabulated WR mass-radius relation of Langer (1989) which 
produces slightly smaller results. 

In the case of a NS secondary, the orbital separation 
does not shrink enough to cast doubt on the results of 
Barnard et al. (2014) and Laycock et ah (2015) in regards 


to the general size of the shadowing absorber and the ex¬ 
tended corona. The binary separation and the orbital veloc¬ 
ities change together, so that the sizes inferred from eclipse 
timing remain relatively unchanged; Table shows a con¬ 
traction of about 20%. The eclipse is still several times larger 
than can be explained by an eclipse of the WR star. Hence, 
the X-ray flux modulation is still required to be dominated 
by absorption and scattering in the dense core of the stellar 
wind, perhaps in the cooler shadowed material (Figure 


Compact object 


Accretion Disk/ Corona 


X-ray Ionized Wind 


Center of Mass -|- 
WRStar 


Shadowing Core of Wind 


Shielded Sec.tor 


The parameter that changes most dramatically is the 
orbital velocity of the WR star vwr^ which decreases by a 
factor of ~10-15 for Mx = 1.4 M© (Table Q. An important 
consequence is that now Vwind » vwr which is a require¬ 
ment for the van Kerkwijk et al.| ( |1996| model to be valid. 
Then the observed RV curve is expected to reflect only a 
small advance in phase due to the orbital and rotational ve¬ 
locities of the WR star that are both very small compared 
to Vwind- The complete picture is more complicated: at in¬ 
creasing distances r from the WR star, the shielded sector 
contains material that was ejected earlier from the point of 
view of an observer corotating in the shadow (Fig[^. So, 
as the wind expands outward, its tangential velocity (fixed 
at the stellar surface value Vrot) is rapidly outpaced by the 
motion of the shadow whose velocity goes as ruj. Thus the 
projected velocity of the emission line will be influenced by 
the phase advance provided by the stellar surface velocity, 
and the trailing which becomes larger at higher altitudes. 
The result will be a line profile that depends on the viewing 
angle and hence on orbital phase. A constant centre of mass 
or 7 velocity is expected in the RV curve, this is not evident 


Silverman V Filippenko (2008) who measured their line 


centroids relative to nebular lines, which probably co-move 
with the binary. In Cyg X-3, 7 has not been measured in 
optical/IR spectra despite that system having a large X-ray 
derived value, a situation Hanson et al. (2000) attribute to 
turbulence. 

Model calculations for IC 10 X-1 by [Tutukov et al.| 


(2013) found that for accepted values of the stellar wind 


parameters in WR stars, the X-ray source should have been 
more luminous by at least one order of magnitude. A simi¬ 


lar conclusion was reached by Barnard et al. (2014), that a 


mere 2% of the Bondi accretion rate is sufficient to power the 
observed X-rays. Both of these calculations assumed large 
BH masses; however, an adopted lower mass (permitted by 
our reinterpretation of the dynamical data) eliminates the 
discrepancy. Specifically, for typical NS parameters {Mx — 
IAMq and Rx — 10 km) and for the measured dynamical 
parameters {Mwr = SBMq, Porbit = 1.45 d, and Vwin d = 
370 km s“^), we find from eq. (4) of Stella et al. (1986) for 
direct wind accretion that Lx = 1.5 x 10^® erg when we 
also adopt the value of M = 2.6 x 10“^ M0yr“^ proposed 
by Tutukov et al. (2013) for the rate of the mass captured 


Figure 1. Schematic diagram of the ICIO X-1 system for M\yji 
= 35 Rq and Mx = 1.4-2.5 Mq with all components drawn 
approximately to scale (See Table 1). X-ray eclipses occur as the 
WR star’s dense wind core region passes through the line of sight 
to the compact object. At mid-eclipse, a “shielded sector” of the 
wind that lies in the shadow of the WR star will be pointing 
directly toward the observer. The vector v^i^d represents the 
velocity of the partially ionized material that gives rise to the 
optical He II emission line. Enhanced absorption and scattering 
by this cooler material contributes to the depth and shape of the 
eclipse minimum. The cross represents the center of mass of the 
binary. The large mass ratio causes the compact object to move in 
a wide orbit at high speed while the WR star moves only slightly 
by comparison. For small Mx, the center of mass lies inside the 
WR envelope. 

3 DISCUSSION 

As a putative 24 M© BH with a 35 M© companion located 
in a low metallicity (Z=0.3 Zq) galaxy, IC 10 X-1 has con¬ 
tributed to the picture of metallicity as determining maxi¬ 
mum stellar BH mass. Dynamically confirmed BHs are few 
in number, so it is striking that the most massive, including 
ULX BHs, and intermediate-mass (IMBH) candidates (e.g. 
HLX-1 in ESO 243-49: Farrell et al.||2014 ) have been found 
in galaxies with extremely low Z. Detailed stellar evolution 


calculations (Belczynski et al. 2010) show that metallicity 
determines the mass-loss rate in the progenitor star’s wind 
and hence the remnant mass. The feasibility of IMBHs grow¬ 
ing in cluster or galaxy cores from stellar BH ‘seeds’ via 
dynamical interactions is critically sensitive to the initial 
stellar BH mass. The above is no longer supported by ICIO 
X-1 if the compact object has a low mass. 

The debate about the nature of the compact objects in 
ULXs (e.g. Gilfanov et al.||2004 Luangtip et al.|[20T4 ) con¬ 
tinues to this day as new observations arise, and the balance 
has tipped recently in favour of stellar-mass BHs and NSs 
Liu et al.||20T3l NGC7 793 P13, [Motch et al.' 


(MlO l ULX-1,^^ ___ _ 

|2014| M-82 X-2, Bachetti et al.|[2014 | rather than IMBHs. 


by the accretor. 


It is actually interesting that the observations do not favour 
the “simpler” resolution of this issue: the main difficulty 
with ULXs is that they should not radiate above the Ed¬ 
dington limit (LEdd = 1.3 X 10^^{Mx/M q) erg s“^), and 
IMBHs with masses Mx^lO^^^M© would remain comfort¬ 
ably below this limit. 

A preliminary attempt to find pulsating ULXs (other 
than M82 X-2) in XMM-Newton archival data ( Doroshenko| 
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et al.||2014 ) found none, but more sensitive surveys will fol¬ 


low. Our attention has turned to ICIO X-1 as a NS candi¬ 
date because of its modest (for a ULX) luminosit y {Lx ^ 
1.2 X 10^® erg s“^) and its phase-shifted RV curve (Laycock 


et al. 2015). In such a case, the Corbet diagram (Corbet 


1986) predicts the NS pulsations to lie in the Pspin<^ s range 


for the observed Por6zt=l-45 d. Notwithstanding the absence 
of any other known WR+NS binary, we suppose that such 
a system will lie between the SG and RLOF systems. 

The lack of pulsations could plausibly be due to lim¬ 
itations of the current datasets. Time resolution in typi¬ 
cal imaging mode for Chandra ACIS is 3.4 s, for XMM- 
Newton EPIC-MOS 2.6 s and EPIC-PN 73 ms. The recent 
discovery of a 1.37 s pulsar in M82 X-2 by |Bachetti et alT] 
(2014) with NuSTAR shows that such objects can hide in 
plain sight if their spin periods lie below the timing thresh¬ 
old. At very short orbital periods (such as Porbit— 1-45 d for 
ICIO X-1 and Porbit^ 2.5 d for M82 X-2), pulsations suf¬ 
fer doppler blurring because of the orbital motion of the NS 
during an observation. The time required for the path length 
variation to accumulate a light-crossing time of Pspinj^ and 
hence to go out of phase is (5t = c PspinI orbit)^ where c is 
the speed of light. This is as little as 250 s for the maximum 
projected velocity in Table of Vorbit =600 km/s (Doppler 
shift z=0.004). At the average uneclipsed XMM-Newton PN 
count-rate of 0.5 c/s, it will be challenging to detect such 
pulsations. 

Stellar winds can wash out high-frequency variability 
when photons reaching the observer scatter off electrons ly¬ 


ing at a wide range of distances from the NS (Kylafis X 
Klimis|1987 ). When the light-travel time across the scatter¬ 


ing region is comparable to the pulse period, the modulation 
is lost. Assuming the scattering core of the WR wind to have 
a radius of ^8 Rq (Table [^, pulse smearing would be signif¬ 
icant for Pspin<^S s. 

A NS with a magnetic field in the usual range for young 
objects (10^^ G) should be a pulsar at IC 10 X-l’s accre¬ 


tion rate. The existence of a slow 7 mHz QPO ( Pasham 
et al. |2013| ) could then be attributed to the interaction of 


the magnetic field with the inner edge of the accretion disk. 
The type of QPO seen in ICIO X-1 is not BH-specific and 
occurs in many NS systems ( Bozzo et al.||2009 ). It has long 
been known that high B-field NSs can switch off their pul¬ 
sations yet continue to emit X-rays. The best known ex¬ 
ample is A0538-66 which has been detected on numerous 


occasions, but only once have pulsations been seen (Skinner 
et al.|1982 |. Magnetospheric accretion that is gated by rota¬ 


tional inhibition may also be able to create unpulsed emis¬ 
sion. This mechanism has previously been invoked to explain 
the persistent low luminosity of some NS+Be systems in qui¬ 
escence and as a way to turn off accretion in SEXTs where 
high mass-transfer rates are continuously present. 


4 CONCLUSIONS 

Based on a finding that the RV curve of the WR X-ray bi¬ 
nary ICIO X-1 is not phase-aligned with the X-ray eclipses, 
we have presented the case that the compact object in this 
system, rather than being a 2O-3OM0 BH, is a NS or a lower- 
mass BH. If the RV curve is a projection of the WR’s wind 
velocity, then there is no longer any dynamical evidence for 


a massive BH in ICIO X-1. In this letter we have focussed on 
application of the van Kerkwijk et al. (1996) model for the 
interaction between the radiation field and the stellar wind 
in HMXBs, and we are keenly aware that the true picture is 
more complex. The wind velocity profile and spatial ioniza¬ 
tion structure remain unconstrained by observation, requir¬ 
ing higher resolution optical/IR spectra as input for model¬ 
ing. The shape of the X-ray modulation is likely related to 
the scattering and absorption properties of the gas in the 
shielded sector which crosses the line of sight at mid-eclipse. 
Hard X-ray observations could probe the effects of scatter¬ 
ing on the eclipse profile, while motivation for phase-resolved 
X-ray line spectra is established by Vilhu et al. (2009). The 
observed X-ray luminosity of X-1 (Wang et al. 2005) lies 
within the range exhibited by both NSs and stellar-mass 
BHs in HMXBs although it is substantially higher than most 
persistent NS+SG systems, and closer to the Be systems at 
the peaks of their giant outbursts. Super-Eddington lumi¬ 
nosity has been seen in A0538-66 and recently M82 X-2, so 
the possibility exists for a NS in IC 10 X-1. The system is 
similar in many regards to other supposed BH-IMXBs which 
now should be examined closely with regard to the possibil¬ 
ity that this phenomenon is more widespread. 

This work was supported in part by NASA grant 
NNX14-AE77G. 
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